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Abstract 

A model is investigated where a monochromatic, spatiaUy homogeneous laser 
field interacts with an electron in a one-dimensional periodic lattice. The 
classical Hamiltonian is presented and the technique of stroboscopic maps is 
used to study the dynamical behavior of the model. The electron motion is 
found to be completely regular only for small field amplitudes, developing a 
larger chaotic region as the amplitude increases. The quantum counterpart of 
the classical Hamiltonian is derived. Exact numerical diagonalizations show 
the existence of universal, random-matrix fluctuations in the electronic energy 
bands dressed by the laser field. A detailed analysis of the classical phase space 
is compatible with the statistical spectral analysis of the quantum model. 
The application of this model to describe transport and optical absorption in 
semiconductor superlattices submitted to intense infrared laser radiation is 
proposed. 
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I. INTRODUCTION 

The search for systems whose dynamical behavior can be fine tuned by one or more 
external parameters has been an important activity in the field of chaos in the last decades 
m] . The kicked rotor has emerged as the paradigm model of a periodically driven system 
whose motion depends strongly on the strength of the external perturbation. The study of 
this particular model allowed us to understand phenomena such as the dynamical localization 
of atoms interacting with microwave radiation in an optical trap . 

For time-independent, Hamiltonian systems, one has also extensively searched for situa- 
tions where the dynamics is a sensitive function of a single parameter. Common situations 
occur when two integrable potentials are coupled to form a nonintegrable system. One case 
of particular interest, both theoretically and experimentally, is the hydrogen atom in the 
presence of a uniform magnetic field 0, since the limiting cases of zero magnetic field and 
zero Coulomb force are exactly solvable. 

In this work we study a model describing the motion of an electron in a one-dimensional 
periodic potential in the presence of an external monochromatic laser field. The most acces- 
sible experimental realization of this model occurs in semiconductor superlattices (multiple 
quantum-well structures) irradiated with intense infrared laser pulses. 

We have considered a one- dimensional time-periodic effective classical Hamiltonian de- 
scribing the interaction of an intense laser field with an electron in a 1-D periodic lattice. We 
show that the effective Hamiltonian is integrable in the limiting case of zero field strength. 
As the field strength becomes larger, the behavior of the classical orbits become increasingly 
chaotic. For strong fields the Hamiltonian is invariant under a scaling transformation involv- 
ing the electron momentum, the average energy of the Hamiltonian, and the field strength. 
Therefore, to understand the whole range of laser intensities it is necessary only to study 
the weak and intermediate field regimes. 

In order to propose a physical realization of this model situation we have studied the 
quantum counterpart of the classical Hamiltonian. Through numerical simulations, we found 
that typical signatures of quantum chaos, such as level repulsion and spectral rigidity, are 
present for intermediate and large values of field strengths. Based on the energy level 
statistics observed, it was possible to identify a weak (Poisson) and a strong (Wigner-Dyson) 
level correlation regime as the laser intensity increases. The statistical spectral analysis 
approach, which provides strong evidence of chaotic and regular behaviors, can be completely 
understood in the light of the classical model. 

This paper is organized as follows. In Sec. we present the classical model and study 
the dynamics of the system as a function of the laser field strength using the technique 
of stroboscopic maps. In Sec. fT^ we derive the quantum model and compare the classical 
results. A discussion on possible experimental observations is given in Sec. |^. Finally, 
Sec. |V] is dedicated to the conclusions. 

II. THE CLASSICAL MODEL 

Consider the Hamiltonian describing an electron in a 1-D periodic potential under the ac- 
tion of an applied homogeneous monochromatic laser field. The laser field can be represented 



by a time-dependent vector potential (in the dipole approximation) minimally coupled to 
the electron. In this limit, the laser field is not affected by the electron motion and the 
vector potential acts as a time-periodic external force. The lattice potential is periodic and 
can be written in terms of a sum over reciprocal vectors Gi = 27ii/d, with d denoting the 
lattice constant and £ = 0, ±1, ±2, .... Explicitly, 
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The term proportional to Aq in H yields a pure time-dependent factor that does not affect 
the equations of motion of the electron. Therefore, this term can be eliminated and the 
classical motion can be obtained from the effective Hamiltonian 

HAp, g) = I^ - 2 £ a, cos {iG,q) + ^^ cos{ut). (2) 
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In order to compare the results of the classical calculations with the quantum analysis 
shown bellow, we first introduce suitable energy and frequency scales into the problem: 
e = h'^G\/2mi and Tiuj = ecuo, such that the dimensionless lattice potential strength becomes 
a^ = ae/e. This naturally leads to a new pair of canonically conjugate variables 

h 
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and a new time scale r = et/h. The resulting equations of motion for P{t) and Q{t) can 
be derived from the rescaled Hamiltonian 

niP, Q) = ^^ - 2 £ a? cos UV2Q) + ^J<^ .osi^or), (3) 

where x = eAoGi/um,c. The total energy associated to the effective Hamiltonian in Eq. (Q) is 
related to the rescaled (dimensionless) energy S{t) by Ees(t) = e £{et/h). The dimensionless 
parameter x characterizes the laser intensity. For numerical estimates, it can be expressed 
in the more appropriate form 
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with / = rirUj'^AQ/STcc and n^ denoting the refraction index of the medium. The quantity 
Ic = m^cu^c^/Svre^ is the Compton intensity, i.e., Ic = 1-37 x 10^^ x A~^ [W/cm^], with the 
laser wavelength A given in /im. 

Since the potential is a periodic function of Q, Q and P are angle and action variables, 
respectively, with P having units of h. Thus, hereafter, we shall drop the h in Eq. (^ and 



understand that P is measured in units of h. We shall also restrict Q to the primitive cell 
interval — 7r/\/2 to 7t/\/2. 

In order to carry out our analysis, we will consider only the first three reciprocal compo- 
nents of the periodic lattice potential in Eq. @. There is nothing special about this choice 
and it just represents a possible truncation of the well-known Kronig-Pennig potential. The 
number of components is taken to be relatively small in order to facilitate the numerical 
simulations; at the same time, it is also sufficiently large to allow for the existence of a 
nontrivial classical dynamics. For similar reasons, in all numerical calculations (classical 
and quantum) in this work we shall adopt a^ = a^ = 0.25 (for —3 < i < 3) and ujq = 0.3. 

Figure |l] displays stroboscopic maps for the classical system defined by Eq. (H) at different 
laser intensities. These maps are generated by plotting the coordinates Q and P at discrete 
times r„ = 27m/uQ, n = 0, 1,2, . . ., producing a sequence of points {Qn,Pn)- Each initial 
condition (Qc-Po) gives rise to a different sequence. If these points fall on a one- dimensional 
curve, the trajectory lies on a cylinder in the extended {Q, P, r) phase space and is said 
to be regular. If, on the other hand, the points cover a two-dimensional area of the QP 
plane, the trajectory is said to be irregular or chaotic. Figure |^ shows contour levels of the 
Hamiltonian for a; = 0. 

Let us now discuss the classical motion arising from Eq. (^ in terms of stroboscopic 
maps for different laser intensities. For x = 0, Figure |l](a), the map is integrable and all 
points lie on 1-D curves. The wells centered at Q = P = and P = 0, Q — ±1.5, trap 
low-energy electrons, whereas those with higher energies, above or bellow these trapping 
islands, move along the lattice from left to right (positive P) or from right to left (negative 
P). These trajectories are clearly related to the conduction properties of the material. In 
Figure |^ the numbers close to the curves indicate the value of S. For x ^ 0, only average 
energies can be defined: above the islands, they are roughly equal to the kinetic term alone, 
since the other terms have an oscillatory behavior. 

For slightly larger intensities, a small amount of chaos begins to permeate the region 
where the separatrices involving the islands used to be. This effect is clearly visible in 
Figure |l|(b), where x = 1. In the new chaotic zone, the electromagnetic field shakes the 
electron, forcing it to move in opposite directions intermittently, inhibiting charge transport 
through the crystal. Figure § shows typical electron trajectories for chaotic (lower curve) and 
regular (upper curve) regimes when x = 1. These trajectories are related to different initial 
conditions. Note that the chaotic trajectory is restricted to a narrower range of Q- values in 
comparison to the regular trajectory, which extends over a much wider Q region. 

For X = 10, Fig .(^), a rather large chaotic zone develops at \P\ < 3.5. Figure |l](d) shows 
that for X = 100, the region close to \P\ = becomes dominated by chaotic trajectories. 
Notice that the maps for a; = 10 and x = 100 are almost identical at large \P\. This is due to 
the existence of a quasi scaling law in the classical motion: For x ^ 1, one can discard the 
periodic potential and change P — > aP and x —>■ ax, rescaling S -^ o?E. As a consequence, 
it is possible to infer that for very strong field intensities the regime will be predominantly 
chaotic, even at (relatively) low energies. 

From the figures and the equations of motion, for any given value of x, it is easy to 
see that P tends to a constant, provided its initial value is sufficiently large. Moreover, at 
X = the average energy also becomes equal to a constant since the average of the coupling 
term becomes negligible. The small P region, on the other hand, is largely affected by the 



electromagnetic field and the lattice potential if x is not too small. 



III. THE QUANTUM MODEL 

The quantum Hamiltonian corresponding to the classical model of the previous section 
can be written in the form 



H = Hk + i/int, (5) 



where 
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(notice that m here is the free electron mass). In these equations we have used Bloch's 
theorem to decouple the Hamiltonian into reduced components, which we denote by the 
label k. We can transform H into a Floquet Hamiltonian Hp by using as a basis the 
eigenstates \i,n) = |G^ + k,n). The Bloch-Floquet states fully incorporate the symmetries 
of the original Hamiltonian: G and n are associate to the discrete space translations of the 
lattice vectors and time translations of the vector potential, respectively. One can then show 
0J^ that the spectrum of H follows from the diagonalization of the Bloch-Floquet matrix 
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where J„ is the Bessel function of order n and ag = a for —3 < i < 3, and zero otherwise. 
The eigenvalues of the Floquet Hamiltonian are quasi-energies but, for simplicity, will be 
referred to as energies. The units and conventions are the same as those defined in Sec. |I|. 
We carry out quantum calculations considering a plane-wave basis set with a finite num- 
ber of Gi : —irnax < ^ < ^maxj which Icads to a (2£max + l)-band model in the absence of 
the laser field. We take ^max = 3 in most calculations presented below, which corresponds 
to a seven-band model. The choice of a finite number of bands to describe a real crystal 
is justified, since highly excited bands (differing from the Fermi level by an amount larger 
that the work function of the material) can never be accessed due to the photoelectric ef- 
fect. In the absence of the laser field, the electronic energy levels of (1^) consist essentially 



of the free-electron parabola folded into the first Brillouin zone, plus gaps opening at the 
crossing points (fc = and ±Gi/2). Inclusion of the laser field produces a dressing effect 
in the bands, which may be described by replica bands translated by nujQ for integer values 
of n. The values adopted for cr° and Uq are such as to strongly mix the replicas forming 
the dressed bands. The interaction ifjnt causes an anticrossing whenever two noninteract- 
ing replica bands cross. Our calculations also involve truncation of the basis set into a 
finite number of n- values: \n\ < rimax- The value of rimax is chosen such as to guarantee 
the completeness of the replica-bands within the energy range of interest, and consequently 
convergence in the calculation of spectral properties (see discussion below) . 

Figure ^ shows the dressed bands spectrum for the energy range l<i5<5ata; = 0, 
X = 1, and X = 10 (as mentioned before, in a real situation the appropriate range of energy 
would be defined by the Fermi energy and the work function of the material). For this 
energy range, n^ax = 40 is sufficient to achieve saturation in the replicas. For a; = no 
level repulsion is observed while for x = 1 some structure in the energy bands is clearly 
identified. For x = 10 the spectrum reveals an intricate level repulsion structure which 
resembles those characteristic of quantum chaotic regimes in other systems |1[]. All the main 
features of the classical stroboscopic map can be inferred from these spectra. For all energies 
we can observe that, as x increases, the anticrossing and enhanced level repulsion among the 
electronic bands leads to narrower minibands. From elementary band theory, narrow bands 
are always associated to poor transport properties. For example, in a standard tight-binding 
scenario, the band width is proportional to the hopping matrix element between atomic 
orbitals (equivalent to Wannier functions in the present plane-waves description) centered 
at neighboring lattice sites. Smaller values of this hopping element indicate a tendency 
towards localization, which, in the corresponding classical motion, appears as trajectories 
confined to a narrower range of Q values (recall similar discussion in Sec. |T| with respect to 
Fig. i. 

These trends are confirmed by a quantitative analysis of the spectral fluctuations. We 
recall [|l| that weak and strong interlevel correlations have been associated with classically 
regular and chaotic behaviors, respectively. The latter is usually described at the quantum 
level by the Gaussian ensembles of random matrix theory [|^. In Figure ^ we present 
the distribution of nearest-neighbor level spacings (NNS) obtained for increasing values 
of the field intensity parameter x, with bands confined to the energy interval 1 to 5. For 
X = [Figure |^(a)] the NNS distribution is very close to an exponential (the Poisson law), 
demonstrating the lack of short-range correlation between levels. The corresponding classical 
momentum varies roughly from -3 to 3 for this energy interval. For x = 1 [Figure |^(b)] a 
crossover regime appears in which the NNS is neither Poissonian nor GOE-like ||^ . However, 
for a; = 10 [Figure ^(c)] the Wigner surmise provides a very accurate fit, as expected from 
the classical analysis. The curves obtained for the least square deviation A3 |0] (not shown) 
lead to a similar interpretation with respect to large-range correlations. 

We now discuss the robustness of our results with respect to the basis set truncation, i.e., 
the influence of the cutoff parameters defining the range of i and n values in the calculations 
presented above. Figure ^ presents the NNS for x = 10. Figure |^(a) is obtained considering 
a plane- waves basis cutoff for £max = 1, while Figure ^b) refers to £max = 5. Note that the 
results in (a) and (b) are essentially the same, and also very similar to Figure ^(c), illustrating 
that the plane-wave basis cutoff is not a relevant parameter. This fact supports that the 

6 



conclusions drawn from our model, where only a finite number of bands is incorporated, 
should be applicable to real experiments. For a given energy range, n^aax must be taken 
sufficiently large to guarantee convergence. Since, in the absence of the laser field, the full 
spectrum associated to larger values of £max covers a wider range of energies, the adopted 
value of rimax increases with £max- Thus n^^x = 60 is required if we consider £max = 5 , 
while for fmax = 3, n^^x = 40 already leads to converged results. The results presented in 
Figure ^ show that the choice of the basis set cutoff parameters must take into account the 
energy range of physical interest. Our main conclusions about the system dynamics inferred 
from the quantum simulations are well understood from the classical analysis. Therefore, 
quantum and classical results are completely compatible. 

IV. EXPERIMENTAL ASPECTS 

The scheme proposed here is suitable to model the physics of high-quality semiconductor 
vertical superlattices in the presence of intense monochromatic light radiation |jlO| . In order 



to make the connection we replace the one-dimensional potential by the superlattice poten- 
tial. It would not be computationally harder to introduce a few more reciprocal components 
to obtain a V{q) in closer resemblance to that of a superlattice. We found, however, that 
increasing the number of components did not change our results qualitatively. Thus, the 
only necessary (and crucial) modification is to replace the free electron mass in Eqs. (^) and 
(0) by the effective mass associated to the material m*. This requires a redefinition of the 
parameter x as well, namely. 




^=h -h-^. (10) 



For a specific application, let us assume a CO2 laser with A = 10.6/im {hu = 117meV) and a 
GaAs-AlrcGai^xAs {x ~ 0.3) heterostructure of lattice constant d = 76A. This choice leads 
to huj = 0.3e. The other parameters take the values m* ~ 0.067m, Ur ~ 3.4, e ~ 390meV, 
a ^ 97meV, and J^ ~ 1.2 x 10^^ W/cm^. We thus find the relation 

9.4x 107 W/cm2' ^ ^ 

Under these conditions, the laser intensity necessary to produce strong coupling, namely, 
X ~ 10, is of the order of lO^'' W/cm^. Such intensity is experimentally accessible and 
does not cause irreversible damage to the heterostructure |]ll|. Nevertheless, one should 



have in mind that not all parameters are fixed by the choice of lattice constant and laser 
frequency. The well width can still be chosen independently (implying in a modification of 
the ratio cr/e). The values chosen above only represent the strong mixing (chaotic) situation 
compatible with our numerical simulations. Of course, decreasing the intensity would make 
the system no longer chaotic. 

We believe that the transition between regular and chaotic behaviors can manifest itself 
in at least two ways. As previously discussed, the onset of chaos is related to a narrowing of 
the electronic bands, causing a decrease in the mobility along the superlattice direction. As 



a result, transport measurements should reveal a maximum of the transverse conductivity 
when the laser intensity is small. Another effect of chaos appears in the strong mixing 
between energy bands. Optical transitions which were initially forbidden when no light is 
shone into the heterostructure should become observable at high laser field intensities. One 
could, in principle, irradiate the superlattice with a pump laser and measure the intensity of 
a large set of these new absorption lines using a second (much weaker) probe laser |0 . The 



histogram of oscillator strengths thus generated should fall into a universal curve predicted 
by random matrix theory |jl2[ when a; ~ 10. 

V. CONCLUSIONS 

We have studied the interaction of a quantized laser field with an electron in a 1-D 
lattice. The classical dynamics of Hamiltonian which describes the problem was analyzed 
through the stroboscopic map technique. We found that a transition from regular to chaotic 
motion can be obtained by increasing the laser field intensity. The main signature of such 
transition would be a large change in the electron mobility, leading to a suppression of charge 
conductance in the irradiated material. An experimental realization of the model is a laser 
interacting with a superlattice. The conditions which must be fulfilled for the laser and the 
lattice parameters are compatible with the experimental possibilities. This implies the use 
of infrared lasers and superlattices of periodicity around lOOA. 
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FIG. 1. Stroboscopic maps of the Hamiltonian in Eq. (^) for (a) x = 0, (b) x = 1, (c) x = 10, 
and (d) x = 100. See the text for a discussion. 
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FIG. 2. Contour plots of the Hamiltonian in Eq. (^) for x = 0. 
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FIG. 3. Classical electron trajectories for x = 1. The upper curve is a regular trajectory for the 
initial conditions Q = 0, P = 3 and propagation time = 50. The lower curve is a chaotic trajectory 
for the initial conditions Q = 0, P = 1.9 and propagation time = 500. 
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FIG. 4. Energy bands for (a) a; = , (b) x = 1, and (c) x = W calculated with parameters 

^max = "J dlid ^max = 40. 
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FIG. 5. Nearest neighbors statistics (NNS) for (a) x = 0, (b) x = 1, and (c) a; = 10 for ^max = 3 
and rimax = 40. The dashed Unes are Poissonian distributions while the soUd hnes correspond to 
GOE distributions of level spacings. 
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FIG. 6. Dependence of the nearest neighbors statistics (NNS) for x = 10 on the number of plane 
waves ^max- The histograms correspond to (a) ^max = 1, ^max = 40, and (b) ^max = 5, nmax = 60, 
i.e., 3 and 11 electronic bands models, respectively. These results, together with Figure |5|(c), 
illustrate the robustness of our results with respect to the basis set truncation parameters. 
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